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I discuss the important role that Lattice QCD plays in testing the 
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1 Introduction 



This year, 2012, marks the beginning of a new era in Physics. A new boson has been 
observed by Atlas and CMS, at a mass of approximately 125 GeV [U [2], that is in 
the range where the Higgs boson, which is the last missing Standard Model particle, 
is expected. Studies on the nature of this particle will be performed aiming at under- 
standing if the Standard Model (SM) is all what we can see in present experiments 
or if New Physics (NP) effects may be revealed. 

In order to search for NP and understand its nature there is a research activity 
that is complementary to the direct production of NP particles, that is Flavor Physics. 
Studies of Flavor Physics look at theoretically clean and SM suppressed processes 
where NP effects may be competitive to the SM contributions and thus visible. Typ- 
ically, the dominant uncertainty in the theoretical predictions of flavor observable 
is brought by the hadronic parameters which enclose the (non-perturbative) long- 
distance QCD contributions. It is therefore crucial to have accurate computations 
of the hadronic parameters. A leading role is played by Lattice QCD as it is a non- 
perturbative approach based on first principles. It consists in simulating QCD itself, 
without any additional parameter, on a discrete space-time and in a finite volume. 

Lattice QCD has recently entered the precision era thanks to the increased compu- 
tational power and the algorithm and action improvements achieved in the last decade. 
The former has led to the so-called unquenched calculations, where the contribution 
of loops of dynamical quarks is taken into account. In the last decade essentially all 
lattice calculations have been performed with either two (up and down) or three (up, 
down and strange) dynamical quarks. Some very recent calculations also include the 
contribution of the dynamical charm quark [3j HI [5] . Thanks to the algorithm and 
action improvements, simulations at light quark masses in the Chiral Perturbation 
Theory (ChPT) regime have become feasible and, very recently, first simulations at 
the physical point have been performed [HI Ej • 

A clear indication of the level of accuracy achieved at present in Flavor Lattice 
QCD calculations is given by the color code introduced by the Flavor Lattice Aver- 
aging Group (FLAG) [8] . The task of FLAG is to review lattice results of interest for 
Flavor Physics and to provide lattice averages, which include lattice results where all 
systematic uncertainties are satisfactorily under control. More in detail, a color tag 
is assigned to the lattice results w.r.t. each systematic uncertainty. Green, orange 
and red colors respectively correspond to the cases of a systematic uncertainty that is 
completely, sufficiently or not enough under control. Lattice results have to have no 
red tags to be included in the average. The first FLAG review j9] provided averages 
for pion and kaon Physics. A second review updating the previous one and including 
also Heavy Flavor hadronic parameters is in progress [10J. A green tag is assigned for 
the continuum extrapolation if the analysis has been performed with at least three 
lattice spacings with at least two values below 0.1 fm. The condition for a green tag 
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for the chiral extrapolation is that the simulated pion masses are lighter than 250 
MeV. The renormalization, where needed, has to be non-perturbative for a green tag. 
Finite volume effects are considered to be completely under control if the product 
^■K,min ' L > 4 or at least three volumes are simulated. The chosen criteria, which 
reflect the state of the art of present lattice results, provide evidence of the high level 
of accuracy achieved in lattice calculations. 

In the following, in order to show the important role of Lattice QCD in Flavor 
Physics, I will discuss an emblematic analysis that relies on several lattice results: 
the determination of the parameters of the Cabibbo-Kobayashi-Maskawa matrix and 
in particular the Unitarity Triangle Analysis performed by the UTfit collaboration. 

2 The Cabibbo-Kobayashi-Maskawa matrix 

One of the main tasks of Flavor Physics is an accurate determination of the parameters 
of the Cabibbo-Kobayashi-Maskawa (CKM) matrix. It represents a crucial test of the 
SM and, moreover, improving the accuracy on the CKM parameters is at the basis of 
NP analyses, where small NP effects are looked for. The CKM matrix, Vckm, is the 
mixing matrix that relates (primed) weak eigenstates to (unprimed) mass eigenstates 
for down-type quarks, as follows 

f d'\ ( d \ ( V ud V us V ub \ ( d \ 

s' = Vckm is \ = \ V cd V cs V cb s . (1) 

V V ) \ b ) \ V td Vu V tb J \ b J 

In the mass eigenstate basis, therefore, the CKM matrix elements appear in weak 
charged currents. Being a 3x3 unitarity matrix, Vckm depends on four independent 
physical parameters: three mixing angles and one phase. In the convenient Wolfen- 
stein parameterization, the CKM matrix is expressed in terms of the four parameters 
A, A, p and rj, as an expansion in the small parameter A, which is the sine of the 
Cabibbo angle (A = sin# c w 0.2). Up to 0(X 5 ), as required by the present level of 
theoretical and experimental accuracy, the Wolfenstein parameterization of the CKM 
matrix reads 

/ 1- \\ 2 - |A 4 A AX^p-nj) \ 

Vckm= \ -A+iA 2 A 5 [l-2(p + ^)] l-iA 2 -|A 4 (l + 4A 2 ) AX 2 .(2) 

V AA 3 [l-(p + ^)(l-iA 2 )] -A\ 2 + ±A(l-2p)\ i -irjAX 4 1 - \A 2 A 4 / 

In the determination of the four CKM parameters Lattice QCD plays a crucial 
role. The parameter A, besides being the CKM expansion parameter, it is particularly 
interesting as it enters the most stringent unitarity condition on the CKM matrix. 
This is, among the nine unitarity conditions V ckm Vckm = 1, the first-row relation 
which reads \V u d\ 2 + \ V US \ 2 + \ V u b\ 2 — 1. In this relation the contribution of \V u b\ 2 can 
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be neglected as it is very small, at the level of and even slightly smaller than present 
uncertainties on \V u d\ 2 and jV^I 2 . The present uncertainty on the first row unitarity 
condition is almost equally distributed between \V UC [\ 2 (4 • 10~ 4 ) and \V US \ 2 (5 • 1CT 4 ). 
The parameter \V u d\ is very precisely determined, at the 0.02% level, from nuclear 
beta decays. It can be alternatively determined, at a similar level of accuracy, from 
the leptonic pion decay, relying on the lattice computation of the pion decay constant. 
The parameter |V^ S | relies on the Lattice results for the kaon decay constant fx or 
for the vector form factor /+(0). The former results allow to determine \ V US \ from the 
experimental measurement of the so-called K12 leptonic decay K — >■ [iv [TTJ, while 
the latter results are required to extract \V US \ from the experimental measurement of 
the so-called K13 semileptonic decay K — > ttIu [121 [13]. F° r a recent review of the fx 
and /+(0) lattice results I refer to the FLAG review [8] and to Gilberto Colangelo's 
plenary talk at the Lattice2012 conference [14J. Here I only quote the FLAG averages 
for \V US \, which combine both the K12 and K13 determinations and are separately 
given for the N f = 2 and N f = 2 + 1 lattice input: \V US \ = 0.2254 ± 0.0009 from 
N f = 2 + 1 and \V US \ = 0.2251 ± 0.0018 from N f = 2. 

As \V US \ is known at present with the impressing precision of 0.5%, small effects 
of the same sub-percent size, like isospin breaking (IB) effects, have now to be in- 
cluded in the determination. So far lattice calculations have been typically performed 
in the limit of exact isospin symmetry, that is with degenerate up and down quark 
masses (m u = rrid) and neglecting electromagnetic effects (Q u = Qd = 0). The 
parametric size of the IB effects is of approximately 1% as they are of 0{a e , m ) or 
0((m d — m u ) I Aqcd) depending on the electromagnetic (Q u ^ Qd) or strong interac- 
tion {m u 7^ m d ) origin. Last year, the strong IB corrections to fx/ fw and to /+(0) 
have been calculated on the Lattice for the first time |15jj The study of ref. [TJ)] is not 
performed removing directly the degeneracy m u = m d , it is instead based on the idea 
of expanding the functional integral in the small parameter {ma — m u )/ Aqcd up to 
first order, with the advantage of computing the (not small) slope in (md — m u )/AQCD- 
By comparing one lattice result, for instance for the kaon mass splitting, to the cor- 
responding experimental value, the quark mass splitting (m d — rn u ) turns out to be 
determined. The results obtained in [15] by adopting this expansion method read 

[m d - m u ] QCD (MS, 2GeV) = 2.35(8) (24) MeV x [ —^ f +J „ , 

L J V ' ; 6.05 x 10 3 MeV 2 

QCD _ M 2 iQCD 

= -0.0039(3)(2) x [Mk ° N \ k +J 2 , (3) 
V A ' 6.05 x 10 3 MeV 2 ' V ; 



Fk/ F n 



where [M^o — Mj c+ ] QCD represents the kaon mass splitting due to strong IB effects. 
In these results the lattice error (the one in the first parenthesis) has been obtained 



*The strong IB effects were previously taken into account in the analysis of [IB] by fitting isospin 
symmetric lattice data through ChPT formulas. 
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with a rather modest statistics, ~ 150 gauge field configurations, and is expected 
to be reduced in the next future. Most of the systematic error (the one in the 
second parenthesis) comes from the ambiguity in the definition of the electromagnetic 
corrections in the experimental input [M^o — M^+]^ CD ', which is also expected to be 
reduced thanks to future lattice computations of the complementary electromagnetic 
effect PHHE]. 

3 The UTA within the Standard Model 

The CKM parameters p and rj are conveniently determined through the so-called 
Unitarity Triangle Analysis (UTA) which consists in constraining sides and 

angles of the triangle defined in the (p, 7/)-plane (p = p(l — A 2 /2) and rj = r/(l — A 2 /2)) 
by the unitarity condition V* b V u d + V* b V C( i + V^ b V t d = which involves the first and 
third rows of the CKM matrix. This triangle has the advantage of having sides of 
similar size and thus of being sensitive to the CP-violating parameter rj. 

Within the UTA several constraints are included, which are provided by the com- 
parison between experimental measurement and theoretical prediction for flavor ob- 
servables that depend on p and rj. The list of the constraints and their present level 
of accuracy is given in table [Tj For some constraints the experimental accuracy is 



Observable 


Accuracy 


£k 


0.5% 


Am d 


1% 


Am d /Am s 


1% 


\Vub/V c b\ 


15% 


Br(B ->■ tu) 


20% 


sin 2/3 


3% 


cos 2/3 


15% 


a 


7% 


7 


14% 


2/3 + 7 


50% 



Table 1: Approximate level of accuracy on the UTA constraints. 

at the level of few percent or even better, so that a significant comparison to the 
theoretical prediction calls for a similarly good control of the hadronic uncertainties. 
The hadronic parameters required in the UTA are the bag-parameter Bk entering the 
theoretical prediction of Ek, the semileptonic form factors /+ and F required for the 
extraction of \V ub \ and \ V cb \ and the combinations of 5( s )-meson decay constants and 
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Figure 1: Results of the UTA within the SM. The contours display the selected 68% 
and 95% probability regions in the (p, rj) -plane. The 95% probability regions selected 
by the single constraints are also shown. 



Observable 


Input value 


SM prediction 


Pull 


e K ■ 10 3 


2.23 ±0.01 


1.96 ±0.20 


1.4 


Am s [ps _1 ] 


17.69 ±0.08 


18.0 ± 1.3 


< 1 


\V cb \- 10 3 


41.0 ± 1.0 


42.3 ±0.9 


< 1 


\v uh \ ■ 10 3 


3.82 ±0.56 


3.62 ±0.14 


< 1 


Br(B tv) ■ 10 4 


1.67 ±0.30 


0.82 ±0.08 


2.7 


sin 2/3 


0.68 ±0.02 


0.81 ±0.05 


2.4 


a 


91° ±6° 


88° ± 4° 


< 1 


7 


76° ± 11° 


68° ± 3° 


< 1 



Table 2: Comparison between input value and SM prediction for the UTA constraints. 
The pull is also shown. 

bag parameters f Bs , fssl Ib, B Bs and B Bs /B B , which enter the theoretical predictions 
of the B-physics observables Am^, Amd/Am s and Br(B — > tv). 

The main results of the UTA [22], performed by the UTfit collaboration assuming 
the validity of the SM, are summarized in fig. [TJ where the curves representing the 
UTA constraints intersect in a single allowed region for (p,rj), proofing that the CKM 
parameters are consistently overconstrained. In other words, the UTA has established 
that the CKM matrix is the dominant source of flavor mixing and CP violation and the 
parameters p and rj turn out to have the values p = 0.139±0.021 and rj = 0.352±0.016. 
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Figure 2: Compatibility plot for ek- The cross, which represents the input value, is 
1.4cr larger than the UTA prediction. 



In table [2] this comparison is shown for all the UTA constraints and the pull (i.e. 
the difference between the input value and the SM UTA prediction divided by the 
uncertainty) is provided as well. For most of the constraints the pull is smaller than 
one, showing that there is a very good compatibility between the input value and the 
UTA prediction. For the three observables Ek, Br(B — > tv) and sin 2/3, instead, there 
is some tension as shown by the pull that is larger than unity and by the compatibility 
plots in figs. |2jj4| 

The theoretical prediction for the CP- violating parameter ek depends on the bag- 
parameter Bk (B denotes the renormalization group invariant B-parameter) which 
encloses the long-distance contribution in K° — K° mixing and for which several 
unquenched results have recently become available. The Nf = 2+1 and Nf = 2 FLAG 
averages read 5^ /=2+1 = 0.738 ± 0.020 and B N K S=2 = 0.729 ± 0.030 The input 
value adopted in the UTA is slightly larger, B K = 0.750 ±0.020, to take into account 
new results [22]- [25] that have appeared after the FLAG review. In particular, the 
result in [23] is characterized by a very safe chiral extrapolation since the simulated 
pion masses are close to the physical value, thanks to the choice of a particularly 
advantageous Lattice QCD action, the 2-step HEX smeared clover-improved Wilson 
action [2"T] . 

The observable where the tension between experimental measurement and UTA 
prediction is the largest is Br{B — > tu), for which the average of the BaBar and Belle 
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BR(B— >xv) 

Figure 3: Compatibility plot for Br{B — » tu). The input value is 2.7a larger than 
the UTA prediction. 




sin2(3 



Figure 4: Compatibility plot for sin 2/3. The input value is 2.4<r smaller than the UTA 
prediction. 
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experimental measurements reads BriB — > Tv) exp = (1.67 ± 0.30) • 10 -4 while the 
UTA prediction, which assumes the SM validity, turns out to be Br(B — > tv)$m — 
(0.82 ± 0.08) • 10~ 4 . In wondering if this 2.7a deviation can be due to NP effects 
the first model that comes to theorists' mind is the simplest 2-Higgs Doublet Model 
(2HDM), that is the 2HDM of type II, where one Higgs boson doublet (H u ) couples to 
up-type quarks and the other Higgs boson doublet (Ha) couples to down-type quarks. 
The observed deviation, in principle, could be easily explained in the 2HDM of type II 
where, in addition to the tree-level SM amplitude with the W-boson exchange, there 
is a tree-level contribution mediated by the charged Higgs. As the Higgs couples more 
strongly to the heavy r lepton than to the lightest muon and electron, the 2HDM of 
type II would seem to provide a natural explanation of the fact that the deviation is 
seen in the r channel only. However, in order to explain the enhancement observed 
in Br(B — > tv), the 2HDM of type II should have a large value of tan P/mjj which is 
instead excluded by other constraints, in particular by the experimental measurement 
of Br(b s'j). At the ICHEP conference an important experimental news has been 
announced by the Belle collaboration [28J. Belle has performed a new analysis, with 
a modified hadronic tag, finding a result for BriB — > tv) that is significantly smaller 
than the previous one and that is compatible with the SM prediction. Thus, the 
present experimental average reads Br(B — > Tv) exp = (0.99 ± 0.25) ■ 10 -4 . Further 
experimental results are certainly looked forward. From the theory side one could 
wonder if the observed enhancement could be due to some underestimated uncertainty 
instead of NP effects. The theoretical prediction for BriB — > tv) is proportional to 
|V^| 2 , which represents the main source of uncertainty in the branching ratio, mainly 
due to the 2.6cr difference between the inclusive and the exclusive determinations of 
this CKM element. The experimental measurements of BriB — > tv) would prefer a 
large value of \V u b\, close to the inclusive determination. However, such a large value 
would not solve rather would worsen the tension in sin 2/3 and, therefore, it does not 
seem to be the solution to the BriB — > tv) puzzle. 

The interest for B decays with a r lepton in the final state has been recently 
stimulated also by the new BaBar (full data) results for the two ratios R(D^) = 
Br(B ->• D^T'Vt) I 'Br(B ->• D^SrVi) [29], which respectively exceed the SM pre- 
dictions by 2.0 and 2.7a, corresponding to a combined discrepancy at the 3. 4a level. 
In two recent papers [301 EI] a more accurate theoretical prediction of the R(D) ratio 
has been provided. The idea of [3U] is to obtain an estimate of R{D) with minimal 
theory input, in particular by using in input the ratio of the vector and scalar form 
factors. In [31], instead, the input value for the scalar form factor is taken from 
unquenched Lattice QCD only [32]. Both papers [30] and [31] slightly reduce the 
discrepancy of the theoretical prediction for R(D) with the experimental measure- 
ment, from 2.0a to 1.8 and 1.7a respectively. The 2HDM of type II, that as in the 
case of Br(B — > tv) in principle could provide an explanation to the enhancements 
in R(D^) in terms of a charged Higgs contribution, would require in this case two 
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different values of to explain the experimental results for R(D) and R(D*). 

More elaborated NP models, instead, could accommodate the enhancements observed 
in Br{B — > tv) and in R(D^). Some of them are 2HDM of type III (where the H u 
and bosons couple to both up- and down-type quarks) with flavor violation in the 
up sector [33] and NP models with right-right vector and right-left scalar currents, like 
some 2HDM, leptoquarks or composite quarks and leptons models (with non-trivial 
flavor structure) [34] . 

4 B-Physics lattice inputs for the UTA 

Lattice results for B-Physics hadronic parameters play a crucial role in the UTA. 
Indeed, the five UTA constraints that rely on Lattice QCD results are Ek, Br(B — > 
tv), Am d , Amd/Am s and \V u b/V c b\, with the last four being B-Physics observables. 

The computation of B-Physics hadronic parameters on the Lattice is complicated 
by the presence of large discretization effects of order (a * irib) up to some power, 
which imply that the physical b-quark mass, being of approximately 4GeV, cannot 
be directly simulated on present lattices (where a -1 < 4GeV). Several methods 
have been investigated and adopted so far, that are either based on an effective the- 
ory approach or consist in simulating heavy quark masses (m/J in the charm region 
(or slightly above) and using some suitable technique to achieve the b-quark region. 
These are the Tor Vergata step-scaling method [35] . which consists in matching sev- 
eral lattice simulations at different volumes and up to the physical b-quark mass on 
the small volume, the ETMC ratio method [3S] based on suitable ratios with exactly 
known static limit, the HPQCD HISQ simulation [37] which reduces discretization 
terms to 0(a s a 2 ml) and C(a 4 m^). Approaches based on effective theory ideas are 
HQET and NRQCD simulations. Within the FermiLab method [38] the key dis- 
cretization errors are removed by tuning three parameters combining perturbation 
theory and experimental inputs. Finally, this year, RBC/UKQCD has formulated 
the non-perturbatively tuned relativistic heavy-quark action [321. ^ * s a variant of 
the FermiLab approach with fully non-perturbative tuning within the By' system. 

In the following of this section I will review the state of the art of the lattice results 
for the B-Physics hadronic parameters entering the UTA. They are the combinations 
of -B( s )-meson decay constants and bag parameters f Bs , fss/ fB, B Bs and B Bs /B B , 
which enter the theoretical predictions of Am^, Am,j/Am s and BriB — > tv) and the 
semileptonic form factors /+ and F required for the extraction of \V u b\ and \V c b\. I 
will also quote the averages, for these hadronic parameters, that are used in the UTA 
by the UTfit collaboration and that represent an update w.r.t. to ref. |40j . 

As far as the decay constants are concerned, it is convenient to consider f Bs , 
which is almost insensitive to che chiral extrapolation as it depends on the light 
quark mass only in the sea, and the ratio f Bs /fB which has the advantage, w.r.t. f B , 
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Figure 5: Collection of recent unquenched results for fs s (left) and fBs/fB (right). 
The red bands represent the average values used as input in the UTA. 



of a partial cancellation of the statistical fluctuations and of the discretization effects. 
Recent lattice results for f Bs and fBs/fB |H]-[11] are collected in fig. [5j The average 
values adopted by the UTfit collaboration have been obtained as simple averages 
of these unquenched (Nf = 2 and Nf = 2 + 1) results, with a conservative error 
that corresponds to the "typical" accuracy of recent calculations. In particular, this 
means that the uncertainty on the fs s average is larger than the error quoted for the 
HPQCD result pE]. The UTA inputs read 

(4) 



f Bs = (233 ± 10)MeV f Ba /f B = 1.20 ± 0.02 
from which it also follows the average for f B given by 

f B = (194 ± 9)MeV . 



(5) 



New accurate analyses are being performed by FNAL/MILC, RBC/UKQCD, ETMC 
and Alpha and have been announced at the Lattice2012 conference [45j-[48j. 

An accurate determination of the B and B s decay constants is of great interest also 
for NP analyses. For instance fs s enters quadratically in the theoretical prediction 
of the branching ratio for the B s — > u, + rare decay, which is highly sensitive to 
NP as it is a loop flavor changing neutral current process and it is theoretically clean 
being a leptonic decay. The experimental bound has been recently strongly improved 
thanks to the LHC measurements. The combination of the the Atlas, CMS and LHCb 
bounds, that is dominated by the LHCb result, provides Br(B s — > < 4.2- 10~ 9 

at 95% CL [49], which is close to the SM prediction. The SM prediction obtained 
from the UTA reads Br(B s — > h + h~)sm = (3.5 ± 0.3) • 10~ 9 and is in agreement 
with the recent result [50], Br(B s — > ^ + [i~)sm = (3.23 ± 0.27) ■ 10~ 9 , that has been 
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obtained in a scheme such that NLO electroweak corrections are likely to be negligible 
and by paying attention to the effect of the soft radiation. 

It is interesting to observe that Lattice QCD can provide useful information also 
for the experimental measurement of some processes like B s — > fi + fi~. Experimen- 
tally, in fact, it is fundamental to know the fragmentation fraction f s /fd that is the 
fraction of b-quarks produced in the hadron collider that hadronize into B s mesons. 
Through factorization, the fragmentation fraction can be related to the ratio of the 
form factors involved in the semileptonic decays B° — > D + £~V and B® — > D^i~V. The 
FNAL /MILC collaboration has computed this ratio on the Lattice (with Nf = 2 + 1 
dynamical fermions and at two lattice spacings) finding (f s /fd) — 0.28 ± 0.04 [32] . 
which represents an improvement w.r.t. to a twenty-years-old QCD sum rule esti- 
mate [5T] . 

The theoretical predictions for B^ — B°^ require, in addition to the decay con- 
stants, the bag-parameters Bb and Bb s - It is convenient, as for the decay constants, 
to take Bbs and the ratio Bb s /Bb in input. For these quantities the UTA lattice 
inputs coincide with the Nf = 2 + 1 HPQCD results |52j 

B Bs = 1.33 ±0.06 B Bs /B B = 1.05 ±0.07, (6) 

as other unquenched results are still preliminary [531 154"] . First unquenched results 
for the bag-parameters of the complete operator basis that describes B^ — B°^ in 
NP models are also looked forward. An analysis by FNAL/MILC is in progress [53] . 

For \V u b\ and \V c b\ there exist two different determinations based on the analysis 
of inclusive or exclusive semileptonic B decays. The inclusive determination is in 
principle less affected by the non-perturbative uncertainties related to the hadronic 
final states. However, as the experimental inclusive measurements require the intro- 
duction of energy cuts, the inclusive determinations of \V u b\ and \V c b\ cannot avoid 
some model dependence in treating long-distance contributions at threshold. This 
is not the case for the exclusive determinations which, instead, rely on theoretically 
clean lattice determinations of the form factors. 

For the exclusive determination of \V u b\ on needs on the experimental side the mea- 
surement of the decay width for B — > ntv and from Lattice QCD the hadronic quan- 
tity T(g 2 > 16GeV 2 )/|Kb| 2 (the larg e-q 2 region is more directly accessible to lattice 
determinations). A collection of the available results for T(q 2 > 16GeV 2 ) / \V u b\ 2 [55J- 
[60] is provided in fig. |6j where the red band corresponds to the average 

iKftlexd = (3.28 ±0.31) -10" 3 , (7) 

which takes into account the (only) two modern unquenched results and older quenched 
results. The comparison to the average quoted by the Heavy Flavor Averaging Group 
(HFAG) [61] for the inclusive determination, |K&|md = (4.41 ± 0.28) ■ 10~ 3 , shows a 
2.6cr discrepancy, indicating that the \V u b\ puzzle is still to be solved. Further lattice 
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Figure 6: Collection of quenched and unquenched results for T(q 2 > 16GeV 2 )/|V^{,| 2 . 
The red band represents the average considered by the UTfit collaboration. 

calculations are certainly desired and are being performed by RBC/UKQCD [62J, 
Alpha [48] and HPQCD [63] or under investigation [64J. The UTfit collaboration 
conservatively combines the two (exclusive and inclusive) values, using as UTA input 
iKtfelinput = (3.82 ± 0.56) ■ 10~ 3 . As we can see from table [2] the UTA prefers a value 
for \V u b\ that is closer to the (lower) exclusive determination. 

The state of the art for V c b presents a similarity to the V u b case. The inclusive 
determination derives from a global fit based on an Operator Product Expansion 
(OPE), in which V c b is fitted together with the b quark mass. The HFAG average [BT] 
reads |V^,|j nc / = (41.9 ± 0.8) • 10 -3 and it is 2.4a larger than the exclusive value 

\V cb \ eX ci = (39.0 ±0.9)- HT 3 . (8) 

We observe that the present accuracy on \ V c b\ is at the 2% level, that is approximately 
five times better than on \V u b\. Both the inclusive and the exclusive determinations, 
in fact, are better under control for \V c b\ than for \V u b\. The reasons are the exper- 
imental cuts at higher energies, where the OPE is more reliable, for the inclusive 
determination, and the fact that the form factors involved in the exclusive determi- 
nation of \V c b\ measure a small deviation from the unity value in the infinite quark 
mass limit. 

Two channels are considered for the exclusive determination, B — > D*iv and B — > 
Div, which respectively require the lattice calculation of the form factors denoted as 
F(l) and G(l). Quenched and unquenched lattice results [65]- [71] for F(l) and G(l) 
are shown in fig. [7j At present the B — >■ D* channel is measured with a better 
accuracy than the B — > D channel, so that the exclusive determination of V c b relies 
on the lattice results for the form factor -F(l). Only one unquenched result [67] 
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Figure 7: Collection of quenched and unquenched results for the form factors F(l) 
(left) and G(l) (right). A star in the legends labels preliminary results. 



exists so far, more recently confirmed by FNAL/MILC itself at Lattice2010 [71] • The 
UTfit collaboration conservatively combines the exclusive and inclusive values, using 
as UTA input |V c f,|j n p Ut = (41.0 ± 1.0) • 10~ 3 . As we can see from table |2j at variance 
with the V u b case, the UTA prefers a value for \V c b\ that is closer to the (higher) 
inclusive determination. 



5 The UTA beyond the Standard Model 

The UTA, besides providing a strong tool for an accurate determination of the CKM 
parameters, can put constraints on possible NP effects. To this purpose the UTfit 
collaboration performs the UTA without assuming the validity of the SM and param- 
eterizing in a model independent way the NP effects that more probably might be 
visible, i.e. the NP contributions in meson-antimeson mixing phenomena (K° — K , 
B° - B°, B° s -B° s ) [72]. Including the new DO data [73J and the recent LHCb mea- 
surement [H] for B° s — B° s mixing, the UTA beyond the SM finds that the NP effects 
in all three systems (K° — K° , B° — B° and B® — B° s ) are constrained to be compat- 
ible with zero [22]- Further measurements are looked forward for the dimuon charge 
asymmetry, in order to confirm or discard the large (non-SM) value of the B° s — B° s 
mixing phase indicated by the DO measurement |75j . 

The NP constraints provided by the UTA analysis beyond the SM can be converted 
into lower bounds on the NP scale. Let us consider, for instance, the K° — K system 
which at present provides the most stringent constraints on NP. In models of physics 
beyond the Standard Model, the effective Hamiltonian that describes the K° — K° 
mixing amplitude involves in general the complete basis of AS = 2 four-fermion 
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operators and it has schematically the form 

^i7 =2 = E^ + E^, (9) 

i=l i=l 

A common choice for the basis is constituted by the operators 

1 = 5^(1-75)^^(1-75)^ , 

2 = ^(1-75)^(1-75)^ , 

3 = s a (l-7s)dV(l-7 5 )d a , (10) 

4 = s Q (l-75)^(l + 75)^ , 

5 = s a (l- 75)^(1 ±7sK , 

where a and /3 are color indices, together with the operators 01,2,3 obtained from 
01,2, 3 with the exchange 75 — > —75. In chirally invariant renormalization schemes, 
the operators Oj have the same matrix elements of the O; and, for this reason, they 
will not be mentioned in what follows. We observe that in the SM only the operator 
0\ appears in the K° — K° amplitude. 

The Wilson coefficients appearing in T-L^ =2 , in an effective theory approach, can 
be parameterized in the form 

Q(A) = ^, i = 2 5, (11) 

where Fi is the (generally complex) relevant NP flavor coupling, Lj is a (loop) factor 
which depends on the interactions that generate Cj(A), and A is the scale of NP, i.e. 
the typical mass of new particles mediating AS = 2 transitions. For a generic strongly 
interacting theory with an unconstrained flavor structure, one expects F,i ~ Lj ~ 1, 
so that the phenomenologically allowed range for each of the Wilson coefficients can 
be immediately translated into a lower bound on A. Specific assumptions on the NP 
flavor structure correspond to special choices of the F^ functions. For example Minimal 
Flavor Violation (MFV) models [TBI [77] correspond to Fx = Fg M and F^i = 0. 

Updated lower bounds on A have been recently obtained in ref. [75], by following 
the same procedure of ref. (79] and using the new unquenched lattice results for B- 
parameters of the complete basis operators calculated by ETMC [78] (with Nf = 2 
and three lattice spacings), which read|T] 

B 2 = 0.54 ±0.03, B 3 = 0.94 ±0.08, B 4 = 0.82 ± 0.05 , B 5 = 0.63 ± 0.07 , (12) 

in the MS scheme defined in ref. [80] at a renormalization scale of 2GeV. The ETMC 
results together with the new RBC/UKQCD results [81] (obtained with Nf = 2 + 1 

^For the definition of the matrix elements in terms of the B-parameters we refer to ref. [75]. 
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Figure 8: Comparison of lower bounds on the NP scale between ref. [78]/ (red) and 
the older analysis J79J (yellow). 

and at one lattice spacing) represent the first unquenched determination for B 2 — B 5 
and turn out to be in agreement. A further computation is being performed by the 
SWME collaboration [82]. The lower bounds on the NP scale as obtained from the 
constraints on the five Wilson coefficients are shown in fig. [8] in a scenario of generic 
flavor structure with tree/strong NP interaction, compared to the older bounds [79] . 
The analysis, which in the considered scenario requires the NP scale to be larger than 
~ 5 • 10 5 TeV, reflects the high sensitivity of Flavor Physics to NP effects. To obtain 
the lower bound on A entailed by loop-mediated contributions, one simply has to 
multiply the quoted bound by a s (A) ~ 0.1 or aw ~ 0.03. 

6 Charm Physics 

In this section I discuss the charm sector of Flavor Physics, which is not involved 
in the UTA and can provide complementary information on the CKM matrix and 
on NP effects. The phenomenon of D° — D mixing, for instance, is sensitive to a 
different sector of NP w.r.t. K and systems since D mesons contains a charm 
quark, which is of up-type and might couple to different NP particles w.r.t strange and 
beauty quarks. From the theory side D° — D mixing has the disadvantage of being 
affected by large long-distance effects, due to the down and strange quarks circulating 
in the loop box diagrams, which dominate over the short-distance contribution. Only 
order of magnitude estimates exist for the long-distance contributions and are at the 
level of the experimental constraints. This prevents the possibility of revealing an 
unambiguous sign of NP. However, barring accidental cancellations between SM and 
NP contributions, significant constraints can still be put on the NP parameter space. 

A very recent analysis of D° — D° mixing, which combines up to date experimental 
measurements, has been performed by the UTfit collaboration [83]. First (prelimi- 
nary) unquenched lattice results for the five i^-parameters entering D° — D° mixing 
in NP models, have been computed with Nf = 2 dynamical flavors by ETMC [M]. 
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Figure 9: Collection of recent unquenched results for f Ds (left) and fo (right). The 
red bands represent the PDG2012 averages. 



They read 

B m = 0.77 ±0.05, B D2 = 0.73 ±0.04, B m = 1.37 ± 0.11 , 

B m = 0.96 ± 0.06 , Bus — 1.22 ± 0.13 , (13) 

in MS at 2GeV. These preliminary values have been used to constrain the parameter 
space of the Minimal Supersymmetric Standard Model (MSSM) with a generic Flavor 
structure [31]. The analysis is performed through the mass insertion approximation, 
that is in the so-called Super CKM basis where squark-gluino couplings are flavor 
diagonal and the squark mass matrix Mj is not. The constraints on the Mj off- 
diagonal elements turn out to be stronger than in the older analysis [HI] by about 
a factor five, mainly thanks to the improved accuracy of the lattice results for the 
^-parameters and the decay constant fr>. 

For the pseudoscalar decay constants fn and frj s several unquenched lattice cal- 
culations have been performed in the last two years. A collection of these accurate 
results is shown in fig. [9j with the red band representing the PDG value [EH] obtained 
by assuming the CKM unitarity. Different lattice determinations turn out to be in 
good agreement and compatible to the unitarity prediction. 

We recall that for f£, s in 2008 there was a deviation at the 3a level between the 
HPQCD lattice result [85] and the PDG experimental average [87] . The HPQCD2007 
result has been superseded by the HPQCD2010 analysis, where the determination of 
the scale has been improved [88], providing fo s = (248 ± 3)MeV that is 2.3cx higher 
than the older result. From the experimental side, new CLEO-c measurements have 
lowered the experimental average by 1.5cr |85j. As a consequence of these two effects 
the Jds puzzle has gone away. For f D , HPQCD has studied the effect of improving 
the scale determination this year [89], finding a results that is in very good agreement 
with the older one [86J and slightly more accurate. First calculations for the D and 
D s decay constants that include the contribution of the charm quark in the sea (i.e. 
with N f = 2 + 1 + 1) are being performed by FNAL/MILC and ETMC. 
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- [20%] 







Table 3: Comparison between the lattice averages for the hadronic parameters enter- 
ing the UTA, quoted by Laurent Lellouch at ICHEP2002 to the values used in input 
by UTfit in the 2012 analysis. Relative uncertainties are shown in square brackets. 

The lattice results for the decay constants /d and Jds represent a fundamental in- 
gredient in the determination of the CKM elements V c d and V cs , respectively, from the 
study of D and D s leptonic decays. At present lattice and experimental uncertainties 
contribute at the same level. An alternative determination of V c d and V cs comes from 
the D — > niu and D — > Ktv semileptonic decays, which require the lattice compu- 
tation of the semileptonic form factors. In this case the present lattice uncertainty 
dominates over the experimental one. The most accurate unquenched results for the 
D — > 7r and D — > K from factors have been obtained by HPQCD [90l |9T] . while a 
new improved analysis is being performed by FNAL/MILC |92j . 

7 A personal outlook to the future 

In the present and next decades there will be a great experimental activity, not only in 
the direct NP searches at LHC, but also in the Flavor sector. Within the quark sector 
the main role in Flavor Physics will be played by LHCb and the SuperB factories. 
The latter experiments aim at improving the accuracy achieved at the B-factories by a 
factor 5 — 10 and, in particular, at testing the CKM matrix at 1% level. They are also 
expected to increase the sensitivity for several channels of interest for NP searches 
by one order of magnitude. Such experimental progress will require the control of 
the theoretical uncertainties, in particular of the lattice uncertainties on the hadronic 
parameters, at the same 1% level. In order to try to understand if such a progress is 
feasible for Lattice QCD I briefly review the progress achieved in lattice calculations 
in the last ten years. In table [3] the lattice averages used in input at present in 
the UTA are compared to the lattice averages quoted by Laurent Lellouch in his 
review talk at ICHEP2002 [93] . The comparison shows that an important progress 
has been achieved in Flavor Lattice QCD in the last ten years, which has typically 
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led to a reduction of the uncertainties by a factor 2 — 5. This has mainly derived 
from the overcome of the quenched approximation, made possible by the increase 
of the available computational power and better algorithms. More recently further 
improvements are being realized, like simulations at the physical point, improved 
control of the discretization effects and the inclusion of the charm quark contribution 
in the sea. I think that we can expect from Flavor Lattice QCD a further significant 
improvement in the next years, toward the 1% accuracy target. 
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